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The monitoring of the gravitational lens Q2237+0305 carried out by the OGLE group 
during 1997-1999 is analyzed. The significant light amplifications in the C and A 
quasar components with maxima in mid- and late 1999, respectively, are interpreted 
as the crossing of microlens caustics by the source. A constraint on the emitting- 
region size R < 10 15 cm has been obtained from the light-curve shape by assuming 
a power- law quasar brightness distribution (r 2 + R 2 )~ p . To estimate the exponent 
p ~ 1.2 requires refining the standard model for the quasar continuum formation in 
an optically thick accretion disk with p — 1.5. 
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1 INTRODUCTION 

Quasar variability under the effect of microlenses depends both on parameters of 
the mass distribution for compact bodies and on the appearance of the emitting 
region. For instance, the larger is the quasar, the smaller are the amplitudes of 
light variations in its images. Thus, it becomes possible to formulate the inverse 
problem of determining the sizes and structure of quasars from their observed 
light curves. The main difficulty in solving this problem is that the distribution 
of the amplification produced by microlenses is not known in advance and ex- 
hibits a fairly complex pattern with many randomly located caustic lines (see, 
e.g., Zakharov 1997 and Zakharov and Sazhin 1998). In general, the specific 
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form of this distribution is not known and can be analyzed only statistically. 
Exceptions are only those portions of the light curve that correspond to caus- 
tic crossing by the quasar. The amplification of a point source during caustic 
crossing obeys a simple law: it remains approximately constant as the caustic 
is approached, then abruptly increases to infinity at the caustic itself, and sub- 
sequently falls off as a; -1 / 2 with increasing distance x from the caustic (Chang 
and Refsdal 1984; Blandford and Narayan 1986). 

Microlenses must show up most clearly in the multiple quasar images pro- 
duced by the gravitational effect of galaxy macrolenses. First, the microlensing 
probability is rather high in such situations. Second, intrinsic quasar variability 
can, in principle, be separated from microlensing variability. Of all the grav- 
itationally lensed quasars, Q2237+0305 (z s = 1.675) is undoubtedly the most 
promising object for microlensing analysis. Because the unique proximity of a 
lensing galaxy (zi = 0.039), microlensing variability in this object must take 
place faster than in other objects by an order of magnitude and with a large 
amplitude. 

The quasar Q2237+0305 has been monitored virtually since its discovery, 
and it actually proved to be the first object in which microlensing variability 
was detected (Irwin et al. 1989; Corrigan et al. 1991). The observations by 
0stcnsen et al. (1996) showed that virtually all four quasar images were more 
or less variable. 

The regularity and quality of the Q2237+0305 monitoring have improved 
markedly when the OGLE group joined it in the last four years (Wozniak et al. 
2000) (see Fig. 1). Measurements are made in the V band approximately once 
a week during the observing season from May through December. The latest 



observations (bttp://www. astro. princcton.cdu/^oglc/ogle2/huchra. html) show 



that image C passed its intensity peak in mid- 1999, while image A peaked in 
late 1999. Interpreting the light-curve maxima as resulting from caustic crossing 
allows the size and structure of the quasar emitting region in the object under 
study to be determined. 

2 AMPLIFICATION OF AN EXTENDED 
SOURCE IN THE CAUSTIC REGION 

When a source crosses a caustic line, an additional pair of images appears (or 
disappears). The total intensity of this pair depends on the distance to the 
caustic as x~ x / 2 . Therefore, the intensity of a point source in the caustic region 
can be represented as (Schneider and Weiss 1987) 

I P (x)=I + 8(x)a x- 1 / 2 . (1) 

Here, Jo is the intensity of all the remaining microimages except the additional 
pair, 9(x) is the Heaviside unit function, and ag is the caustic strength. The 
amplification of an extended source with a brightness distribution P(r) can be 
calculated by ordinary summation over the set of infinitely small sources with 



2 



' .0 [ 



0.8 - 



0.6 - 

I ;a 

hi 

0.4 - 



0.2 - 



600 



800 




1000 1200 

JD 24500004 



1400 



1600 



1800 



Figure 1: V light variations of the four images of the quasar Q2237+0305 during 
1997-2000 as observed by the OGLE group 



individual amplification factors: 
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(2) 



To describe the brightness distribution in the source, we use a power-law model 

P(r) = (l + ^)- p (p>0), (3) 

which is determined by the source radius R and by the rate of brightness decline 
P- 

Let the source center be at a distance D from the caustic line. In the 
normalized coordinates £ = x/R and r\ = y/R and using normalized distance 
d = D I R, we obtain 



I = I + a R- 1 ' 2 J{d), 



where the function J{d) is 
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or, taking into account the symmetry in 77 and the total intensity 
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Figure 2: Function J(d) for several exponents p: curve 1,2, and 3 for p = 1.5, 
2.0, and 3.0, respectively 



we obtain: 

J(d) = 2 P — dd d V cH 1 + (t-d) 2 + r 1 2 T P - (7) 
K Jo Jo 

Taking the internal integral over r] yields 

Ad) = ^ jf Q.f - ^ r 1/a [e 2 - 2£d + + 1] ^ p , (8) 
where B is the beta function. The subsequent integration over £ yields 

r(p-§)r(2p-f)„, , j2 ^_„ „ „ 3 i ^ , d 
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Here, T and 2-Fi(o, 6; c; z) are the gamma function and the Gauss hypergeomet- 
ric function, respectively. Figure 2 shows the function J(d) for several values of 
p. We clearly see from the figure that the sharpness of the jump in amplitude 
increases with increasing source brightness concentration toward the center dur- 
ing caustic crossing, tending to an infinite point-source limit for very large p. 

For some particular p values, for example, for p = 3/2, 2 and 3, the hyper- 
geometric function can be expressed as a combination of elementary and other 



4 



special functions. Thus, we have for p = 3/2 



2(l + rf 2 )((l + d 2 ) 1 / 2 -d) 



1-1/2 



(10) 



3 VARIABILITY ANALYSIS 



We see from Fig. 1 that the intensities of all four quasar images have varied 
during the last four years. The largest variations were observed in image C, 
which passed its maximum in mid- 1999, and in image A, which reached its 
maximum in late 1999. The approach to interpreting the variability of these 
two images is the same. Let us assume that, in both cases, the source crosses 
the caustic; the brightness distribution must follow the law (Q). The general 
form of the curve depends on five parameters: 

1. Contribution Iq from the remaining microimages, which is assumed to be 
approximately constant; 

2. Caustic strength do; 

3. The time it takes for the source to cross its radius At of the source, which 
is proportional to the source size R; 

4. The time to of caustic crossing by the source center; 

5. Exponent p in the brightness distribution (^|). 

Estimating the five parameters reduces to minimizing the sum of the squared 
differences between model and observed light curves 



The summation is performed over all N points of observations with weights that 
are inversely proportional to the squares of the observational errors <7j. The best 
set of parameters is sought by the Marquardt method (see, e.g., Chapter 15.5 
in the book by Press et al. 1992). 

3.1 Image A 

The model parameters estimated from the light curve of image A with 182 
data points during the entire observing period 1997-1999 are given in the Table 
1. The combination ao/At is more convenient to calculate than the caustic 
strength clq. The listed formal accuracies of the parameter estimates should be 
considered only as their lower limits. 

Figure 3 shows the model light curve together with measured values. The 
1997 observations are poorly fitted by a single curve, implying that approaching 
the caustic did not show up in the first year. Such a behavior is characteristic 
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Image 


Io 


a /At 


At 


to 


P 


A 1997-1999 


0.44±0.01 


0.73±0.07 


91±4 


1462±1 


1.24±0.03 


C 1998-1998 


0.22±0.01 


0.95±0.15 


-29±2 


1390±1 


1.10±0.02 



Tabic 1: Best-fit model parameters 



0.90 p 



0.80 - 



0.70 - 



0.60 



0.50 




500 



1000 1500 
JD 2450000+ 



2000 



Figure 3: Variability of image A during 1997-2000. Observational data, model, 
and forecast 



of numerical microlensing models for Q2237+0305, in which one amplification 
event is often superimposed on another to form complex light curves. Exclud- 
ing the 1997 data from our analysis, while significantly improving the total \ 2 
residual, affects the parameter estimates only slightly. 

When the caustic is crossed in image A, an additional pair of images emerges. 
A characteristic feature of this direction of motion is a steep rise in light am- 
plification followed by a gentler decline. The dashed line in Fig. 3 represents 
the expected behavior of image A in the immediate future. We assume that 
the brightness will continue to decline and (if no additional causes of amplifi- 
cation arise) will reach the original 1997 level in two to three years. The latest 
observations for 2000 are added for comparison. 

3.2 Image C 

Attempts to fit the light curve of image C over the entire observed period have 
failed. At the same time, excluding 80 data points for 1997 from our analysis 
results in quite reasonable estimates, which are given in the table and shown in 
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Figure 4: Light curve of image C during 1997-2000. Observations, model (dot- 
ted line), and forecast 

Fig. 4. Interestingly, the best solution corresponds to caustic crossing in the 
negative direction, with the pair of images disappearing. The local minimum of 
X 2 corresponding to the motion in the positive direction is several times greater 
than the absolute minimum reached during the motion in the negative direction. 

The inability to fit the entire observed period by a single curve becomes 
more understandable in light of the recent results by Wyithe et al. 2000. These 
authors argue that the intensity variations in image C have been composite in 
pattern during the last four years, and, what is quite possible, another caustic 
crossing was overlooked during the observing seasons in 1997 and 1998. 

The dotted and dashed lines indicate the best extrapolation computed in 
the model with single caustic crossing in 1997 and our forecast until the end of 
2000, respectively. No appreciable intensity variations in image C are expected 
in the immediate future. 

4 THE SIZE AND STRUCTURE OF 
THE QUASAR EMITTING REGION 

For the quasar radius estimates At in time units to be transformed to linear 
sizes R, we must know the apparent quasar velocity; only the velocity compo- 
nent v± is perpendicular to the caustic line is important. Of course, the exact 
velocity is unknown. Nevertheless, a statistical analysis of the time derivatives 
of brightness variations by Wyithe et al. 1999 yielded an upper limit of v < 500 
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km/s. The perpendicular velocity component can only be lower than this value. 
The most probable value of v± was found to be 300 km/s. Given the angular 
distances from the observer to the quasar D os and to the lens D„i in the model 
of a flat Universe with H = 60 km/s Mpc~ x , the most probable source radius 
is estimated to be 

R = v ± At^ ~ 2.0 ■ 10 13 ( — ^ ; ) ( cm. (12) 
D oi \300km/sJ \lday J K ' 

Using the crossing time At a = 90 days for image A gives the most probable 
quasar radius R = 1.8 • 10 15 cm, while substituting the crossing time for image 
C, Ate = 30 days, reduces this estimate to R = 6 • 10 14 cm. At the same time, 
the velocity constraint v < 500 km/s together with the crossing time Ate give 
an upper limit on the source radius, R < 10 15 cm. 

Another parameter that characterizes the mass distribution is the exponent 
p. It follows from the table that its computed value lies between 1.1 and 1.25. 
It is interesting to note that different models for the quasar structure can lead 
to different dependences of the emissivity on distance from the center. For 
instance, a power-law dependence with p < 0.5 follows from the model of an 
optically thin accretion disk (Manmoto et al. 1997). At the same time, the 
standard model of an optically thick accretion disk yields an r~ 3 dependence 
(Shakura and Sunyaev 1973), which changes to (r 2 + i? 2 )~ 3 / 2 for a finite radius. 
Our estimate p <~ 1.2 favors the standard model, but more complex accretion- 
disk models should be analyzed to achieve better agreement. 

5 DISCUSSION 

Let us consider the legitimacy of some assumptions made here. The main as- 
sumption is associated with the hypothesis of caustic crossing. Numerical calcu- 
lations show that there are two effects capable of causing a significant increase in 
amplification. Apart from fold-caustic crossing, the source can also pass near the 
caustic cusp. However, events of the second type for the images of Q2237+0305 
are several times less probable than those of the first type (Wambsganss et 
al. 1992; Lewis and Irwin 1996) and, in general, are more symmetric. These 
two properties can serve as a statistical justification for using the hypothesis of 
caustic crossing. 

A power-law model with a finite core radius is used to calculate the quasar 
size and structure. Three models of a caustic-crossing source are encountered 
in the literature: a homogeneous disk, a Gaussian source (Schneider and Weiss 
1987), and a (1 - r 2 /i? 2 ) 1/2 distribution (Schneider and Wagoner 1987). All of 
them are completely determined by their radii. The (1 + r 2 / R 2 )~ p model differs 
radically in that it is a two-parameter model. Estimating the rate of brightness 
decline p allows us to choose between different models for the quasar structure. 

Allowance for the spatial orientation of the emitting region appreciably com- 
plicates the analysis. Two additional parameters associated with the orientation- 
angle components of the initial circular source appear. The first and second 
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parameters determine, respectively, the apparent-ellipse eccentricity and the 
angle of motion with respect to the perpendicular to the caustic (in this case, 
allowance for the direction of motion does not reduce to a simple substitution 
of the transverse velocity). Besides, accretion-disk rotation can cause an addi- 
tional asymmetry of the emitting region in the spectral range considered due to 
the Doppler effect. 

In view of many influencing factors, the possibility of reconstructing, at 
least in principle, the one-dimensional quasar profile (along the x-axis) as the 
solution of an integral equation (Grieger et al. 1991; Mineshige and Yonehara 
1999; Agol and Krolik 1999) seems of great interest. However, this is a separate, 
independent problem, which is yet to be applied to an actual monitoring. 

The simple amplification behavior <~ x~ x / 2 for a point source located in 
the caustic region is possible only when several conditions are satisfied. The 
radius of caustic curvature must be considerably larger than the source size (see 
Fluke and Webster 1999 for curvature allowance). The caustic must be isolated 
lest the source cover several caustics simultaneously. In addition, it is implied 
that there is no large-scale time gradient in amplification variations. Thus, for 
example, introducing a constant slope as an additional free parameter allows a 
model curve of image C to be easily fitted to the entire observing period 1997- 
1999. However, since there were no such large gradients throughout the entire 
15-year-long monitoring history of Q2237+0305 and since their physical origin 
is not completely understood, we have to abandon the additional parameter. 

Finally, the problem can be further complicated by intrinsic quasar vari- 
ability. In general, the latter is rather difficult to separate from microlensing 
variability. However, given that the delays between the images in our cases 
do not exceed several days, intrinsic quasar variability must be repeated in 
all quasar images (with individual amplification factors) virtually simultane- 
ously. The fact that the brightness variations in all four components are not 
synchronous provides circumstantial evidence that intrinsic quasar variability is 
negligible in this source. 

Despite the possible complicating factors listed above, we have every reason 
to believe that the simple model used here is capable of faithfully reproducing 
the observational data, and that its implications deserve a careful analysis. 

After this paper was mainly complete, Wyithe ct al. 2000 independently 
published a preprint where they also interpreted the OGLE-group observations 
of Q2237+0305. These authors analyzed the light curves by using statistical 
methods based on conditional probability distributions. 

Wyithe et al. 2000 focused mainly on computing the probability of oc- 
currence of brightness bursts with observed parameters and on estimating the 
possibility of subsequent bursts. The conclusion that there is an additional over- 
looked brightness burst in image C associated with caustic crossing between the 
observing seasons of 1997 and 1998 seems of considerable interest. Such an 
event can account for the difficulties of modeling the light curve for image C in 
terms of the model with a single caustic crossing. Wiythe et al. 2000 also expect 
another caustic crossing in image C 500 days after the 1999 summer maximum 
(with a large uncertainty of ~ 100 — 2000 days, though). 
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The amplification in image A in late 1999 is interpreted as caustic crossing 
in the negative direction, just as we did here. However, the peak of image C in 
mid-1999 is assumed to be caused by the passage of the source near the caustic 
cusp. The choice between the two interpretations of image C variability could 
be made by analyzing color variations of the source, which are much larger 
during caustic crossing than during cusp passage. A color analysis could be 
performed by invoking additional data of the simultaneous monitorings at the 
Apache Point and Maidanak Observatories (Bliokh et al. 1999) through various 
filters during the same observing period. 
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before their publication. 
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